Quantum efficiency of a microwave photon detector based on a current-biased 

Josephson junction 
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We analyze the quantum efficiency of a microwave photon detector based on a current-biased 
Josephson junction. We consider the Jaynes-Cummings Hamiltonian to describe coupling between 
the photon field and the junction. We then take into account coupling of the junction and the 
resonator to the environment. We solve the equation of motion of the density matrix of the resonator- 
junction system to compute the quantum efficiency of the detector as a function of detection time, 
bias current, and energy relaxation time. Our results indicate that junctions with modest coherence 
properties can provide efficient detection of single microwave photons, with quantum efficiency in 
excess of 80%. 
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Quantum optical photodetectors have contributed sig- 
nificantly to the development of quantum optics and 
atomic physics [l| and now play an essential role in op- 
tical quantum information applications such as quantum 
key distribution [Hi)]. Recently, a novel paradigm for 
the study of radiation-matter interaction based on cir- 
cuit quantum electrodynamics (cQED) has emerged as a 
test-bed for quantum optics in mesoscopic systems Q . 
In addition, cQED architectures are also an attractive 
candidate for scalable quantum computing and transmis- 
sion of quantum information 0, H| . Following the orig- 
inal proposal, a variety of cQED architectures demon- 
strating strong coupling between single photons and su- 
perconducting i nteg rated circuits have been realized ex- 
perimentally |9|, |lOj. This work has paved the way for 
the development of a superconducting microwave photon 
detector [ll], EH with possible applications to quantum 
information processing and communication [l3| . 

The microwave photon detector is based on the 
current-biased Josephson junction (JJ): the JJ is biased 
so that absorption of a single microwave photon induces 
a transition to the voltage state, producing a large and 
easily detected classical signal [12[. While these detec- 
tors are straightforward to operate and show potential 
for scalability, performance is degraded by spurious dark 
counts due to quantum tunneling events in the absence of 
an absorbed photon; moreover, energy relaxation within 
the detector results in photon loss and leads to a reduc- 
tion in measurement efficiency. 

In this work, we theoretically determine the quantum 
efficiency of a microwave photon detector based on a 
current-biased J J. Previous analysis [H, 14 of this sys- 
tem was focused on a continuous flux of photons incident 
on the detector. Here we study the probability of pho- 
ton detection by a JJ coupled to a microwave mode that 
is loaded with a fixed number of photons. We solve the 
equation of motion for the density matrix of a resonator 
coupled to a JJ to obtain detector efficiency for differ- 
ent values of operation time, current bias, and relaxation 



FIG. 1: (a) Schematic circuit diagram of JJ-based microwave 
photon detector coupled to a resonator, (b) Potential energy 
landscape of the detector when bias current is close to the 
critical current of the JJ . The junction is initialized in the 
\g) state and upon absorbing an incident photon transitions 
to the \e) state, which rapidly tunnels to the continuum. 



time of the junction. Our results indicate a JJ with decay 
time around 10 ns can detect a single photon in the mi- 
crowave mode with efficiency greater than 80% for read- 
ily achievable circuit parameters. We also find that the 
detector efficiency increases significantly with increase in 
energy relaxation time T± of the junction, suggesting that 
a highly efficient single microwave photon detector is at- 
tainable for a moderate improvement in junction quality. 

The circuit diagram of the JJ and resonator system is 
shown in Fig. QJa). The J J is biased with a current I 
close to its critical value Iq. The junction Hamiltonian 
can be written in terms of the charge operator Q and the 
operator 5 of the phase difference across the J J [l5| : 



U(S) 



2tt 



cos 6 6 

Iq 



(1) 

Here C is the junction capacitance, <l>o = 27r?i/2e is the 
magnetic flux quantum, and h and e are Planck's con- 
stant and the electron charge, respectively. For I < I 0l 
the potential energy landscape U(S) takes on a "tilted 
washboard" shape, with a few discrete energy levels in 
shallow minima separated from the continuum by a bar- 
rier, see Fig.QJb). We truncate the junction Hamiltonian 
to the ground \g) and first excited states |e) and obtain 
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FIG. 2: (Color online) Switching probability Pi (Pq) us. pho- 
ton detection time for initial states \l,g) (\0,g)). Parame- 
ters used in this plot are: junction Ti = 10 ns, barrier height 
AU/fkjp = 2, vacuum Rabi frequency Q./2n — 100 MHz, and 
cavity decay time 1/k = 1 /is. Here the detuning between 
the resonator and the junction A = 0. The solid blue curve 
is the quantum efficiency r)i = Pi — Po (see the text). The 
maximum quantum efficiency of the detector is 50% for the 
optimal detection time of 40 ns (the optimal point is marked 
by a filled blue circle in the plot). In the inset, we plot the 
quantum efficiency vs. detuning A/fi for the optimal detec- 
tion time td = 40 ns obtained at A = 0. The bandwidth of 
the detector is around 3.4 tt. 



tion is modeled by the Jaynes-Cummings (JC) Hamilto- 
nian 1161 : 
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where f2 is the vacuum Rabi frequency and <r+ (<r_ ) is the 
junction creation (annihilation) operator. 

Our goal is to describe the time evolution of the den- 
sity matrix p(t) for the resonator and the junction. We 
represent the density matrix in the basis where 
n measures the number of photons in the resonator, and 
j G {g, e} represents the ground or excited state of the 
junction. We describe coupling of the system to its envi- 
ronment in terms of phenomenological decay rates. The 
coupling of the resonator with the environment produces 
a resonator decay rate k — u r /Q, where Q is the quality 
factor of the resonator. Similarly, finite dissipation in the 
J J induces decay of the excited state with rate 7 ~ 1 /RC 
(decay time T\ — 1/7), where R represents the effective 
damping in the junction. The time evolution of the den- 
sity matrix p(t) is governed by the following equation: 



dp{t) 
dt 
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the following Hamiltonian for the J J: 



H 



jj 



(2) 



where H e = [e)(e| is the projection operator to the ex- 
cited state and uj eg = (e e — e s )/S. 

The Hamiltonian of a single mode microwave resonator 
with resonant frequency uj r is given by 



H r = hoj r (a f a + 1/2) , 



(3) 



where (a) is the creation (annihilation) operator of the 
resonator. The coupling of the resonator with the junc- 



where C-y[p(t)] and C K [p(t)] are superoperators that cap- 
ture damping in the junction and the resonator at low 
temperatures T <C fko eg , fvjj r (l7j : 

As[p(*)] = K ^a/5a f - ]^{a ] a, p]^ , (6a) 

£7 = 7 - \{v+°-,p} \ ■ (6b) 



To account for switching of the ground and the first 
excited states of the junction to the voltage state 
introduce the tunneling superoperator £t[/5(0] 0, llSf 
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where r e . 9 are the tunneling rates from the ground, \g), 
and first excited, |e) states of the junction. 

The system is originally prepared in a pure state 
Pn(0) = \n,g){n,g\ with n photons in the cavity and 
the junction in the ground state \g). We numerically 



solve the above equation for time evolution of the den- 
sity matrix to compute occupation probabilities of the 
resonator and junction states. The probability that the 
JJ has switched to the voltage state at time t is given by 

P n {t) = 1-Tr [p n (t)}. (8) 
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FIG. 3: (Color online)(a) Quantum efficiency rji vs. photon 
detection time for A U /ftw p = 2 and for various decay times of 
the junction: T\ = 10 ns (solid blue), Ti = 20ns (solid black), 
Ti = 50 ns (solid green), Ti = 100 ns (dotted black), and Ti = 
500ns (dotted red). For decay time of 500 ns, the maximum 
quantum efficiency is 94%. (b) Quantum efficiency 771 vs. 
photon detection time for junction Ti = 10 ns for various bias 
points of the junction: AU/fkJ p = 2 (solid blue), 1.9 (solid 
black), 1.8 (solid red) and 1.7 (dashed black). The maximum 
quantum efficiency is 84% for AU/fiUp — 1.7. The remaining 
parameters are as in Fig. [2] 
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FIG. 4: (Color online) Efficiency r} n to detect a photon vs. 
detection time for n = 1 (solid blue), 2 (dashed blue), and 
3 (dash-dotted blue) photons in the resonator. The rest of 
the parameters are as in Fig. [2] For three photons in the 
resonator, the efficiency to detect a photon is 85%. 



We consider the following set of parameters for the 
JJ-resonator system, typical of those realized in ex- 
periments characteristic frequency of the junction 



/2tt = 4.8 GHz, junction decay rate 7 = 10 s , res- 



onator decay rate k = 10 s , and vacuum Rabi fre- 
quency 0/27T = 100 MHz. The tunneling rate of the 
first excited state of the junction is approximated by: 
7.3 x lO^" 1 " 



T P « 500 r 



21] 



In Fig. [2J we plot switching probabilities P\ (t) , Pq (t) of 



the junction for initial states \l,g) (solid red) and |0,<?) 
(solid black), respectively. For this part of the simula- 
tion, we consider the parameter AU/hujp = 2, where oj p 
is the plasma frequency of the JJ, and the detuning be- 
tween the resonator and the junction is set to zero, i.e., 
A = u r — Lj eg = 0. The switching probability Pi in 
Fig. [2] features steps whose periodic occurrence is a man- 
ifestation of Rabi oscillations of the junction with period 
7r/f2 = 5 ns. This result is consistent with the picture 
that switching of the junction halts momentarily when 
the junction returns to the ground state in course of the 
Rabi oscillations. 

Next, we discuss the presence of the wide plateau of 
Pi in Fig. [21 The occurrence of this plateau can be un- 
derstood from the fact that switching of the junction is 
briefly frozen when the junction relaxes to the ground 
state due to dissipation. The JJ then switches to the 
voltage state after time ~ l/T g , the characteristic time 
scale for switching of the junction in the case of zero 
photons. The height of this plateau can be estimated as 
r e /(r e + 7 ) « 0.5, which agrees with the numerical result 
in Fig. H 

In order to determine the quantum efficiency of the de- 
tector, we must properly treat dark counts due to quan- 
tum tunneling from the |0, g) state in the absence of pho- 
ton absorption. The quantum efficiency 771 of the detector 
is defined as the difference between the switching proba- 
bilities for an initial state with one photon P\(t), and for 
an initial state with no photons Po{t): r\\ = Pi(t)—Po(t). 
The quantum efficiency is shown in Fig. [2] by the solid 
blue curve. For our choice of parameters T e ~ 7, the 
detector has maximum efficiency of about 50% for the 
optimal detection time td around 40 ns. 

Next, we examine the bandwidth of the Joscphson mi- 
crowave photon detector. Here, we vary the frequency Lo r 
of the resonator and compute the quantum efficiency of 
the detector for the optimal detection time td obtained 
at zero detuning A = 0. The detector bandwidth is then 
given by the detuning at which the quantum efficiency 
of the detector reduces to half the efficiency obtained at 
zero detuning. For a dissipation-free junction, the de- 
tector bandwidth is approximately given by the vacuum 
Rabi frequency. However, in the presence of dissipation 
and tunneling, the first excited state of the junction is 
broadened by ~ 7 + r e . This broadening of the energy 
level roughly accounts for the increased bandwidth of the 
detector. We find that the bandwidth is a factor of 3.4 
larger than the vacuum Rabi frequency for the param- 
eters of the detector chosen for Fig. [21 as shown in the 
inset of Fig. [21 see inset. 

Next, we analyze the effect of dissipation and bias point 
on the efficiency of the detector. In Fig.(3ja), we plot the 
quantum efficiency of the detector for different values of 
the junction relaxation time T\ from 10 ns to 500 ns, 
keeping all other parameters the same as in Fig. [21 We 
find that the quantum efficiency increases with increasing 
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junction relaxation time T\ and reaches 94% for T\ = 
500 ns. 

The change in bias current / of the junction modifies 
the ratio of barrier height AU to the junction plasma 
frequency u p . Taking different values of this ratio, we 
compute the efficiency of the detector at fixed relaxation 
time T\ = 10 ns; the results are shown in Fig.[3](b). Upon 
decreasing the ratio AU/hu> p , the potential well becomes 
shallower, leading to enhanced tunneling of the first ex- 
cited state. This in turn increases the efficiency of the 
detector. Our simulation results indicate that a signif- 
icant improvement in the detector efficiency is achieved 
when the tunneling rate exceeds the dissipation rate of 
the junction. We find that for bias point AU/Hu) p = 1.7, 
the efficiency of the detector is about 84% for detection 
time around 9 ns. 

Finally, we analyze the efficiency of the detector when 
the resonator is loaded with n > 1 photons. Generalizing 
the previous case of a single photon, the efficiency to de- 
tect a single microwave photon in resonator loaded with 
n photons is given by r\ n = P n (t) — Po{t)- In Fig. 01 we 
plot the efficiency at fixed bias point AU /hu p = 2 and 
T\ = 10 ns for different numbers of photons in the res- 
onator: n = 1 (solid blue), 2 (dashed blue), and 3 (dash- 
dotted blue). We find that detection efficiency increases 
upon increasing the number of photons in the resonator 
and reaches 85% for three photons in the resonator. This 
result is consistent with previous studies 12|, ll4j that re- 
ported higher efficiency of the detector, for the same pa- 
rameters as above, when a continuous flux of photons was 
incident on the detector. For the case of a single photon 
in the resonator, the detector returns to the ground state 
after the photon is absorbed by the environment, and no 
further excitation of the junction is possible. However, 
when multiple' photons are present in the- resonator, other 
photons are available to induce excitation if the junction 
relaxes, thereby increasing the probability of photon de- 
tection. 

In conclusion, we have presented a model to determine 
the quantum efficiency of a microwave photon detector 
based on a current-biased J J. We have demonstrated that 
the efficiency to detect a single photon loaded in the res- 
onator has maximal value r e /(r e + 7). We have also 
determined that the bandwidth of the detector is charac- 
terized by the sum of the vacuum Rabi frequency and the 
broadening of the first excited state of the JJ due to tun- 
neling and relaxation processes. Our simulations indicate 
that for currently used J J photon detectors, the quantum 
efficiency is about 50% for the bias point AU/hw p = 2 
and about 85% for AU/hu> p = 1.7. We have finally in- 
vestigated the probability to detect a photon in the case 
of a multiphoton initial resonator state and have found 
that the detection efficiency quickly approaches 100% as 
the initial number of photons increases, consistent with 
previous studies 12|, Il4j of a continuous flux of photons 
incident on the detector. 
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